Eulerian Rotations of Slater Determinants on a 3-D Cartesian Lattice 
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We discuss three practical methods for performing Eulerian rotations of Slater determinants in a 
three-dimensional Cartesian geometry. In addition to the straightforward application of the active 
form of the quantum mechanical rotation operator, we introduce two methods using a passive 
position-space rotation followed by an active spin-space rotation, one after variation and the other 
before variation. For increased numerical accuracy, the computations are done using basis-spline 
functions. 
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I. INTRODUCTION 

The description of a nucleus in terms of a self- 
consistent mean-field is one of the most commonly used 
approaches in nuclear many-body calculations [3, 0]- 
However, this approach breaks many of the fundamental 
symmetries present in the nucleus, thus making it desir- 
able to restore some of these symmetries and consider 
the configuration mixing of symmetry-restored mean- 
field states This allows the study of correlations in 
addition to those present in the standard mean-field re- 
sults (e.g. Pauli correlations). The correlated states are 
usually generated by symmetry restoration via the pro- 
jection of angular momentum, parityjparticle number (in 
case of pairing), and others 0, H, H, 0, Hi • 

Since the most common choice for a many-body wave- 
function appropriate for mean-field calculations is a 
Slater determinant, many of the projections mentioned 
above involve some application of the rotation operator 
onto a Slater determinant. For example, angular momen- 
tum projection from a Slater determinant representing a 
deformed nucleus involves integration over three Euler 
angles of all possible rotations of the Slater determinant. 
The standard method used in doing such rotations in an 
unrestricted three-dimensional geometry is not only very 
CPU intensive but can also suffer from accuracy prob- 
lems [9]. In the past, these problems were circumvented 
by only dealing with axially symmetric nuclei, which re- 
sults in reducing the integration to being over a single Eu- 
ler angle. Even in this situation a simple finite-difference 
discretization of the derivative operators used in rotation 
may not provide sufficient accuracy. 

In this manuscript we discuss three practical methods 
for performing Eulerian rotations of Slater determinants 
in a three-dimensional Cartesian geometry. In addition 
to the straightforward application of the active form of 
the quantum mechanical rotation operator, we introduce 
two methods using a passive position-space rotation fol- 
lowed by an active spin-space rotation, one after varia- 
tion and the other before variation. The rotation before 
variation approach can also be utilized to generate exact 
initial states for time-dependent Hartree-Fock (TDHF) 
calculations involving deformed nuclei at various orien- 
tations. For increased numerical accuracy, the computa- 



tions are done using basis-spline functions. Basis-splines 
have been effectively used in HF [n|, TDHF [ll|, and 
HFB calculations El El 



II. ACTIVE AND PASSIVE ROTATIONS 

In the quantum theory of rotations, it is important to 
distinguish carefully between active and passive transfor- 
mations 



14l | ■ In an active rotation the wavefunction of 



a quantum particle is rotated by an angle 9, and the co- 
ordinate system is left unchanged. In a passive rotation 
the coordinate system is rotated by an angle —9, and the 
wavefunction of a quantum particle is left unchanged. 

In the following we consider an active rotation of a 
spinor wavefunction (spin 1/2) by an angle 9 around an 
axis described by the unit vector n. The rotation opera- 
tor R has the structure 0, El 



R(9,n) 



exp 



-i— n • J 

n 



(1) 



where J is the total angular momentum operator, con- 
sisting of orbital and spin angular momenta 



-ih(r x V) 



(2) 



and cr denotes the vector formed by the three Pauli ma- 
trices. Instead of performing a rotation by an angle 
9 around the axis given by the vector n, one usually 
parametrizes the rotation in terms of three Euler angles 
(a,/?, 7). First, one makes a rotation through an angle 
7 about the original z-axis, then a rotation through an 
angle (3 about the original y-axis, and finally a rotation 
through an angle a about the original z-axis. Thus the 
operator for an active rotation has the structure 

R{Sl) = R(a, 13, 7) = e -^Jjn e -i0j v /h e -nJ z /h ^ 

We insert Eq. @ into Eq. © and observe that rotations 
in position space, generated by L, and rotations in spin 
space, generated by S, commute. This leads to 



R(n) = R s (n) Rr(Q) 



(4) 
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with the spin space rotation operator 

R„(Sl) = e -^/2 e -</9^/2 e -i7^/2 t ( 5 ) 

and the position space rotation operator 

R r (Q) = e -i<xL*/1' e -i0Ly/h e -i 1 L„/h _ ( 6 ) 

We use the general expression for the spin rotation oper- 
ator [13] 

JZ s (0, n) = e -<«nWa = cos ^ i s _ j s i n _ n . a j (7) 

where I s denotes the 2x2 unit matrix in spin space. Us- 
ing this expression to calculate the operators in Eq. l[5]). 
we obtain the desired form of the spin rotation operator 
in terms of the Euler angles 




,e + ssm&e 2 e + 2cos^-e + 2 



We now apply the rotation operator, Eq. to a given 
spinor wave function tp to obtain the rotated spinor wave 
function ip' 



The expression given in Eq. ([9]) is the most common 
approach for rotating Slater determinants by rotating 
single-particle spinors. However, one can formulate an 
alternate approach where the single-particle spinors are 
unchanged but the coordinate system is rotated instead. 
This corresponds to a passive rotation in position space, 
and it can be expressed in the form [l6| 



Rr(Q) 



^i(r) 
V>2(r) 



^i(r') 
V<2(r') 



(10) 



with 



r' = R-^fl) r 



(11) 



The operator denotes the inverse of the operator 

for an active rotation in position space or, equivalently, 
the operator for a passive rotation of the coordinate sys- 
tem. This operator is given by [ijj 



R-\n) 



COS7 sin7 0> 
— sm.7 COS7 
1, 

cosa sina ' 
— sina cosa 
1, 



'cos/3 —sin/3' 

1 
sin/3 cos/3 



(12) 



(9) or, explicitly [13] 



R;\n) 



/ cos7COS/3cosa — shrysina cos7COS/3sina + sin7Cosa — cos7sin/3> 
-sin7COs/3cosa — cos7sina — sin7COS/3sina -I- cos7COsa sin7sin/3 



sin/3cosa 



sin/3sina 



3/3 / 



(13) 



In summary, we can compute the rotated spinor wave 
function ip via the operations 




where the spin rotation matrix is given in Eq. ||SJ) and 
the rotated position vector is given in terms of Eq. (fTTj) 
and Eq. 

Generalization of the above transformations to a many- 
body Slater determinant is straightforward due to the 
one-body nature of the involved operators. It can be 
shown that the passive transformation for a Slater deter- 



minant $ for a A-nucleon system is simply [1 



JZ r (n)tf(ri,ra s ... ) r A ) = 

*(i? I 7 1 (r!)r 1 ,i?- 1 (r!)r 2 ,...,i?- 1 (r!)rA) , (15) 



followed by a similar transformation in spin-space. We 
also note here that if the calculations are done in a basis 
comprised of analytic functions the implementation of 
such rotations will become straightforward. 
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III. APPLICATION AND RESULTS 

In the sections below, we discuss three methods that 
can be used to perform an Eulerian rotation of a Slater 
determinant and provide relative accuracies for the meth- 
ods. The first two methods deal with rotations per- 
formed after the variational solution of the Hartree-Fock 
(HF) equations. The other method involves the initial- 
ization of the HF calculations with a rotated initial guess 
and the subsequent minimization. A detailed description 
of our three-dimensional, unrestricted HF code is given 
in Refs [13, EH- For the effective interaction we have 
used the Skyrme SLy4 force [isj |. Since all the calcula- 
tions are done with no geometrical symmetry restrictions, 
the single-particle spinors carry only state indices that 
count the state number and isospin quantum number. 
For coordinate-space lattice solutions of the mean-field 
equations, a very accurate representation of derivative 
operators and functions is necessary for performing ac- 
curate rotations. This may be understood by viewing 
the rotation as an interpolation onto a rotated coordi- 
nate system, whereby the mapping of the original lattice 
points to the new points results in the interpolation of 
the single-particle states, which requires high accuracy 
to preserve the integrity of the system. Basis-spline in- 
terpolation [lj| provides one such accurate approach, and 
we use it in our HF calculations. For all of our compu- 
tations, we have tried to maintain the numerical details 
consistent with performing realistic calculations. For test 
purposes we used a deformed 20 Ne nucleus on a 30 3 lat- 
tice with a box size of (—10, +10) fm in each direction. 



A. Rotation After Variation 

1. Active Rotation via Eq. [3| 

This approach is the commonly- used method for imple- 
menting the action of the rotation operator onto a Slater 
determinant. The rotation is done after the minimum- 
energy solution to the HF equations is achieved. For 
Euler angles (a, /3, 7), we apply Eq. ^ to each of the 
single-particle states directly by breaking up the action 
of the exponential operators into small angular steps 
(Aa, A/3, A7), where each step is carried out by the ex- 
pansion of the exponential as a Taylor series. This yields, 
for example, for a rotation by angle a about z, the re- 
currence relation 191 



> = 



JV 

E 

71=0 



(-iAaJ z y 



(16) 



The choice of N = 4 seems adequate for angular steps 
on the order of Aa = 27r/360. The operations given in 
Eq. (|16p require the action of various derivative operators 
onto the single-particle states. For example, for the case 



TABLE I: Rotation for Euler angle (3 = 45° as a func- 
tion of basis-spline order M for the rotation method of Sec- 
tion MI A II The exact values for the tabulated quanti- 
ties are: E B = -157.2688 MeV, r rms = 2.9236 fm, e n = 
-13.2430 MeV, and e p = -9.3527 MeV. 





M=3 


M=5 


M=7 


M=9 


Eb 


-165.6985 


-157.3751 


-157.2654 


-157.2640 


Trms 


2.9558 


2.9258 


2.9237 


2.9236 




-12.8422 


-13.2192 


-13.2406 


-13.2419 


e p 


-9.0123 


-9.3320 


-9.3497 


-9.3508 



above, J z has the differential form 



Jz = ih(yd x - xd y )I + -a z 



(17) 



Similar expressions are for the rotation of (3 about y 
and that of 7 about z. In the Basis-Spline Colloca- 
tion Method (BSCM) fl9| . the derivative operators for 
each coordinate become matrices represented on the lat- 
tice and operate on the discretized single-particle states. 
Thus we are able to compute the above expansion far 
more accurately than can be done with finite-difference 
discretization. 

In testing this method, we have found that the maxi- 
mum error occurs for 45° rotations for each of the Euler 
angles and is worst for angle (3. The reason is that this 
orientation requires the most significant interpolation ac- 
curacy relative to the original code frame (although we 
are not doing interpolation here, the discretization of a 
differential equation can be viewed as an interpolation 
problem) . In Table Q] we show the results for the 45° (3 
rotation as a function of the basis-spline order M. Tabu- 
lated are the total HF binding energy, the r.m.s. radius, 
and the single-particle energies for the least-bound neu- 
tron and proton states. These results are obtained by 
recalculating the HF Hamiltonian using the rotated set 
of single-particle states. The results show the cumula- 
tive error for the (3 rotation as is done from 0° to 45°. 
The spline order M = 3 would be equivalent to using a 
three-point formula for discretizing the derivative oper- 
ators. As we can see, the results for higher-order spline 
functions give reasonable accuracy for this rotation. In 



TABLE II: Rotation with basis-spline order M = 9 us- 
ing the rotation method of Section IIII A II Here we are 
varying a, (3, and 7. The exact values for the tabulated 
quantities are: E B = -157.2688 MeV, r rms = 2.9236 fm, 



e„ = -13.2430 MeV, and e v 



-9.3527 MeV. 



(45,45,45) (90,90,90) (135,135,135) (180,180,180) 

E B -157.2607 -157.2582 -157.2532 -157.2586 

r rm3 2.9236 2.9235 2.9235 2.9235 

e n -13.2414 -13.2401 -13.2401 -13.2403 

e„ -9.3498 -9.3474 -9.3455 -9.3477 



Table [H] we tabulate the cumulative errors for rotating 
in all three Euler angles for basis-spline order M = 9. 
Again, while not perfect, the errors are reasonable for this 
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set of numerical parameters. The difference of the direct 
application method discussed in this section relative to 
the ones discussed below is that the errors accumulate 
due to the recursive nature of the computations. 



requirement is comparable to or better than the direct 
active method mentioned earlier, depending on the num- 
ber of points used in the discretization of the Eulcr angles 
(more points is favorable here). 



2. Passive Rotation by Spline Interpolation 

In the BSCM approach, each component of the single- 
particle spinor is expanded in terms of basis-spline func- 
tions of order M 

1>(x, y,z)=J2 C ijk Bf I (x)Bf(y)Bjf(z) , (18) 

ijk 

where the Cafe's are expansion coefficients in terms of the 
basis-spline knots, the B's are basis-spline functions, and 
we have omitted the single-particle state indices for nota- 
tional simplicity. In the basis-spline collocation method, 
x, y, and z are then discretized on a collocation lattice 
such that each component of the single-particle spinors 
attains the form 



ijk 



Bf{x x )B™{y,)B™{ Zv 



3 



which, for a given order M, may be expressed in matrix 
notation 



E 



Cij kB{\ Bj ^B^n 



(19) 



where A, and v are the indices for the three- 
dimensional Cartesian collocation lattice. This expan- 
sion leads to the corresponding lattice representation of 
the Hartree-Fock equations [13, Ell- After we obtain 
the collocation-lattice solutions of the Hartree-Fock equa- 
tions in the unrotated system, we can compute the ex- 
pansion coefficients in Eq. (fl"9|) as [l^] 



C ijk =Y,B Xi BMB» k Tj }Xllv , 



(20) 



where B denotes the inverse of Bi\ and so on. Conse- 
quently, the passive rotation operation given by Eq. (jTTj) 
simply reduces to an interpolation of the single-particle 
states at the rotated coordinate points 

4{x', y', z') = ]T C ijk B™{x')B¥{y')B*f{z') . (21) 

ijk 

We then rotate the resulting spinors in spin space via 
Eq. ([8]). The expansion coefficients Cijk need only be 
calculated once for each single-particle spinor component 
and can be used for as many (a, (3, 7) combinations as 
are needed. The calculation of the expansion coefficients 
is the most CPU-expensive part of this procedure and 
has the same dimensions as the wavefunction array in 
the code. However, because we need to calculate the 
expansion coefficients only once, the effective CPU time 



TABLE III: Rotation with basis-spline order M = 9 us- 
ing the rotation method of Section IIII A 21 Here we are 
varying a, /3, and 7. The exact values for the tabulated 
quantities are: E B = -157.2688 MeV, r rms = 2.9236 fm, 
e„ = -13.2430 MeV, and e v = -9.3527 MeV. 

(45,45,45) (90,90,90) (135,135,135) (180,180,180) 



E B 



-157.2587 
2.9235 
-13.2410 
-9.3488 



-157.2688 
2.9236 
-13.2428 
-9.3525 



-157.2595 
2.9235 
-13.2411 
-9.3488 



-157.2688 
2.9236 
-13.2428 
-9.3525 



In Table |HI| we again tabulate the cumulative errors for 
rotating in all three Euler angles for basis-spline order 
M = 9. What is interesting for this case is that the 
errors do not accumulate for larger angles because we 
are not reaching those angles recursively but by direct 
interpolation. As a matter of fact, every 90° rotation 
for each angle reproduces the exact solutions due to the 
symmetry of the mesh in the three Cartesian directions. 
In this sense, if an integral over the three Euler angles 
were considered (like in angular momentum projection) 
the cumulative error is expected to be much smaller using 
this passive transformation. 



B. Rotation Before Variation 

In a fully-unrestricted, three-dimensional calculation 
for a deformed nucleus, all orientations of the nucleus are 
equivalent. In practical calculations the final orientation 
of the nucleus depends on its initial orientation. In other 
words, the filling of the levels for the initial guess for 
the HF single-particle states determines the final orien- 
tation. Normally this guess results in a deformed nucleus 
whose symmetry axis (say an axially-symmctric nucleus) 
is oriented along one of the principal axes of the code 
frame. However, by rotating the initial guess via the pas- 
sive transformation given by Eq. (fTTj) and a subsequent 
rotation in spin-space, we can obtain solutions to the HF 
equations that are oriented exactly at those angles with 
respect to the unrotated solution. This method of rota- 
tion before variation yields the most accurate result of 
all three methods. The only discrepancies again seem to 
arise for 45° rotations, most likely because all the deriva- 
tive operators are generated along the principal axes of 
the code frame and may be slightly inefficient when act- 
ing on a nucleus which is maximally-rotated with respect 
to the principal axes. While this method may be more 
CPU-intensive if solutions are needed at many different 
values of the Euler angles, it is very suitable for the ini- 
tialization of time-dependent Hartree-Fock (TDHF) cal- 
culations for deformed nuclei [13, HH . Since the initial 
HF states are usually analytic, the computation of the 
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states on the rotated coordinates is trivial and the con- 
verged HF state is an eigenstate of the HF Hamiltonian 
without any approximations, which is the most suitable 
initial state for TDHF calculations. In Table ITVl we show 
the same cumulative errors shown for the previous two 
methods. As we can see, this method is considerably 
more accurate for intermediate angles and is exact for 
90° rotations for each angle. 



TABLE IV: Rotation with basis-spline order M — 9 using 
the rotation method of Section IIIIBI Here we are varying 
a, /3, and 7. The exact values for the tabulated quanti- 
ties are: E B = -157.2688 MeV, r rms = 2.9236 fm, e„ = 
-13.2430 MeV, and e p = -9.3527 MeV. 

(45,45,45) (90,90,90) (135,135,135) (180,180,180) 

E B -157.2679 -157.2688 -157.2679 -157.2688 

r rm s 2.9234 2.9236 2.9234 2.9236 

e„ -13.2430 -13.2430 -13.2430 -13.2430 

e v -9.3523 -9.3527 -9.3523 -9.3527 



stricted calculations do not preserve the quantum num- 
bers associated with the nuclear system under study, it is 
often desirable to project various good quantum numbers 
from these determinants. These projections frequently 
involve the rotation of the determinental wavefunction 
as the main computational ingredient. In this work, we 
have presented a comparative study of three such meth- 
ods, the first being the standard active rotation employed 
in most calculations. The other two methods involve a 
passive rotation in coordinate space followed by an ac- 
tive rotation in spin space. We have discussed the ad- 
vantages and disadvantages of each method and have 
demonstrated that the accurate implementation of three- 
dimensional rotations do require increased accuracy in 
representing the differential operators on the lattice. The 
basis-spline expansion provides one such approach, which 
has been demonstrated by the accuracy obtained for all 
three methods. Unfortunately, we do not have a finite- 
difference code to compare our results with directly. Our 
method for doing rotations before variation is very useful 
for initializing TDHF calculations for deformed nuclei at 
different orientations. 



IV. CONCLUSIONS 



With the advances in computational power, it is in- 
creasingly possible to perform nuclear structure and reac- 
tion calculations without any spatial symmetry assump- 
tions. One subset of such calculations involves the mean- 
field methods with the principal wavefunctions being a 
single or multiple Slater determinants. As such unre- 
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